


**Figure 1


ttest dis if wave==1 & inlist(q1random, 1, 2), by(q1control)
matrix mat1=(r(mu_1)-r(mu_2), r(mu_1)-r(mu_2)-1.96*r(se), r(mu_1)-r(mu_2)+1.96*r(se))*100

ttest dis if wave==2 & inlist(q1random, 1, 2), by(q1control)
matrix mat2=(r(mu_1)-r(mu_2), r(mu_1)-r(mu_2)-1.96*r(se), r(mu_1)-r(mu_2)+1.96*r(se))*100

ttest dis if inlist(q1random, 1, 2), by(q1control)
matrix mat3=(r(mu_1)-r(mu_2), r(mu_1)-r(mu_2)-1.96*r(se), r(mu_1)-r(mu_2)+1.96*r(se))*100

matrix mat4=(0.16, 0.1, 0.22)*100


***

ttest neu if wave==1 & inlist(q1random, 1, 2), by(q1control)
matrix mat5=(r(mu_1)-r(mu_2), r(mu_1)-r(mu_2)-1.96*r(se), r(mu_1)-r(mu_2)+1.96*r(se))*100

ttest neu if wave==2 & inlist(q1random, 1, 2), by(q1control)
matrix mat6=(r(mu_1)-r(mu_2), r(mu_1)-r(mu_2)-1.96*r(se), r(mu_1)-r(mu_2)+1.96*r(se))*100

ttest neu if inlist(q1random, 1, 2), by(q1control)
matrix mat7=(r(mu_1)-r(mu_2), r(mu_1)-r(mu_2)-1.96*r(se), r(mu_1)-r(mu_2)+1.96*r(se))*100

matrix mat8=(-0.04, -0.09, 0.02)*100


**

ttest app if wave==1 & inlist(q1random, 1, 2), by(q1control)
matrix mat9=(r(mu_1)-r(mu_2), r(mu_1)-r(mu_2)-1.96*r(se), r(mu_1)-r(mu_2)+1.96*r(se))*100

ttest app if wave==2 & inlist(q1random, 1, 2), by(q1control)
matrix mat10=(r(mu_1)-r(mu_2), r(mu_1)-r(mu_2)-1.96*r(se), r(mu_1)-r(mu_2)+1.96*r(se))*100

ttest app if inlist(q1random, 1, 2), by(q1control)
matrix mat11=(r(mu_1)-r(mu_2), r(mu_1)-r(mu_2)-1.96*r(se), r(mu_1)-r(mu_2)+1.96*r(se))*100

matrix mat12=(-0.12, -0.17, -0.08)*100

**

matrix tomz = mat4\mat8\mat12
matrix rownames tomz = a b c 

matrix wave1 = mat1\mat5\mat9
matrix rownames wave1 = a b c 

matrix wave2 = mat2\mat6\mat10
matrix rownames wave2 = a b c 

matrix wave12 = mat3\mat7\mat11
matrix rownames wave12 = a b c 

**

coefplot (matrix(tomz[.,1]), ci("tomz[.,2] tomz[.,3]") label(Tomz (2007)) ms(Oh)) ///
 (matrix(wave1[.,1]), ci("wave1[.,2] wave1[.,3]") label(Wave 1) ms(Dh))  ///
 (matrix(wave2[.,1]), ci("wave2[.,2] wave2[.,3]") label(Wave 2) ms(Th))  ///
 (matrix(wave12[.,1]), ci("wave12[.,2] wave12[.,3]") label(Waves 1 and 2) ms(Sh)),  ///
 coeflabel(a = "Disapproval" b = "Neutral"  c = "Approval") ///
mlabel format(%9.2g) mlabposition(1) xline(0, lp(dash)) graphregion(color(white)) ///
ciopts(recast(rcap)) legend(row(1) pos(6))


***Figure 2: by level of esclation

ttest dis if wave==1 & inlist(q1random, 1, 2), by(q1control)
matrix mat1=(r(mu_1)-r(mu_2), r(mu_1)-r(mu_2)-1.96*r(se), r(mu_1)-r(mu_2)+1.96*r(se))*100

ttest dis if wave==2  & inlist(q1random, 1, 2), by(q1control)
matrix mat2=(r(mu_1)-r(mu_2), r(mu_1)-r(mu_2)-1.96*r(se), r(mu_1)-r(mu_2)+1.96*r(se))*100

ttest dis if inlist(q1random, 1, 2), by(q1control)
matrix mat3=(r(mu_1)-r(mu_2), r(mu_1)-r(mu_2)-1.96*r(se), r(mu_1)-r(mu_2)+1.96*r(se))*100

matrix mat4=(16, 10, 22)

ttest dis if wave==1 & inlist(q1random, 1, 3), by(q1control)
matrix mat5=(r(mu_1)-r(mu_2), r(mu_1)-r(mu_2)-1.96*r(se), r(mu_1)-r(mu_2)+1.96*r(se))*100

ttest dis if wave==2  & inlist(q1random, 1, 3), by(q1control)
matrix mat6=(r(mu_1)-r(mu_2), r(mu_1)-r(mu_2)-1.96*r(se), r(mu_1)-r(mu_2)+1.96*r(se))*100

ttest dis if inlist(q1random, 1, 3), by(q1control)
matrix mat7=(r(mu_1)-r(mu_2), r(mu_1)-r(mu_2)-1.96*r(se), r(mu_1)-r(mu_2)+1.96*r(se))*100

matrix mat8=(16, 10, 22)

ttest dis if wave==1 & inlist(q1random, 1, 4), by(q1control)
matrix mat9=(r(mu_1)-r(mu_2), r(mu_1)-r(mu_2)-1.96*r(se), r(mu_1)-r(mu_2)+1.96*r(se))*100

ttest dis if wave==2  & inlist(q1random, 1, 4), by(q1control)
matrix mat10=(r(mu_1)-r(mu_2), r(mu_1)-r(mu_2)-1.96*r(se), r(mu_1)-r(mu_2)+1.96*r(se))*100

ttest dis if inlist(q1random, 1, 4), by(q1control)
matrix mat11=(r(mu_1)-r(mu_2), r(mu_1)-r(mu_2)-1.96*r(se), r(mu_1)-r(mu_2)+1.96*r(se))*100

matrix mat12=(23, 16, 29)


ttest dis if wave==1 & inlist(q1random, 1, 5), by(q1control)
matrix mat13=(r(mu_1)-r(mu_2), r(mu_1)-r(mu_2)-1.96*r(se), r(mu_1)-r(mu_2)+1.96*r(se))*100

ttest dis if wave==2  & inlist(q1random, 1, 5), by(q1control)
matrix mat14=(r(mu_1)-r(mu_2), r(mu_1)-r(mu_2)-1.96*r(se), r(mu_1)-r(mu_2)+1.96*r(se))*100

ttest dis if inlist(q1random, 1, 5), by(q1control)
matrix mat15=(r(mu_1)-r(mu_2), r(mu_1)-r(mu_2)-1.96*r(se), r(mu_1)-r(mu_2)+1.96*r(se))*100

matrix mat16=(32, 26, 39)

**

matrix tomz = mat4\mat8\mat12\mat16
matrix rownames tomz = a b c d

matrix wave1 = mat1\mat5\mat9\mat13
matrix rownames wave1 = a b c d

matrix wave2 = mat2\mat6\mat10\mat14
matrix rownames wave2 = a b c d

matrix wave12 = mat3\mat7\mat11\mat15
matrix rownames wave12 = a b c d

**

coefplot (matrix(tomz[.,1]), ci("tomz[.,2] tomz[.,3]") label(Tomz (2007)) ms(Oh) ) ///
 (matrix(wave1[.,1]), ci("wave1[.,2] wave1[.,3]") label(Wave 1) ms(Dh))  ///
 (matrix(wave2[.,1]), ci("wave2[.,2] wave2[.,3]") label(Wave 2) ms(Th))  ///
 (matrix(wave12[.,1]), ci("wave12[.,2] wave12[.,3]") label(Waves 1 and 2) ms(Sh)),  ///
 coeflabel(a = "Threat of Force" b = "Display of Force"  ///
 c = "Use of Force without Casualty" d = "Use of Force with Casualty") ///
mlabel format(%9.2g) mlabposition(1) xline(0, lp(dash)) graphregion(color(white)) ///
ciopts(recast(rcap)) legend(row(1) pos(6))

**Appendix A

**Table A1
sum age age30 education college male ccp international coastal west central if wave==1

**Table A2
sum age_group age30 college male ccp international coastal west central income socialstatus if wave==2

**Appendix B

**Table B1
tab q1random, gen(scenario)

logit scenario1 age education male ccp international coastal west central regime goal cost interest if wave==1
logit scenario2 age education male ccp international coastal west central regime goal cost interest if wave==1
logit scenario3 age education male ccp international coastal west central regime goal cost interest if wave==1
logit scenario4 age education male ccp international coastal west central regime goal cost interest if wave==1
logit scenario5 age education male ccp international coastal west central regime goal cost interest if wave==1

**Table B2

logit scenario1 age_group college male ccp international coastal west central income socialstatus regime goal cost interest if wave==2
logit scenario2 age_group college male ccp international coastal west central income socialstatus regime goal cost interest if wave==2
logit scenario3 age_group college male ccp international coastal west central income socialstatus regime goal cost interest if wave==2
logit scenario4 age_group college male ccp international coastal west central income socialstatus regime goal cost interest if wave==2
logit scenario5 age_group college male ccp international coastal west central income socialstatus regime goal cost interest if wave==2

***


**Appendix C

gen ratings = 8 - q1response
replace ratings = 6-q1rating if wave==2

**Table C1:

reg ratings scenario2-scenario5 if wave==1
outreg2 using wave1, excel replace 

reg ratings scenario2-scenario5 regime goal cost interest if wave==1   
outreg2 using wave1, excel append 

reg ratings scenario2-scenario5 regime goal cost interest age education male ccp international coastal west if wave==1  
outreg2 using wave1, excel append 


**Table C2:


reg ratings scenario2-scenario5 if wave==2
outreg2 using wave2, excel replace 

reg ratings scenario2-scenario5 regime goal cost interest if wave==2   
outreg2 using wave2, excel append 

reg ratings scenario2-scenario5 regime goal cost interest age_group college male ccp international coastal west income socialstatus if wave==2
outreg2 using wave2, excel append 


**Appendix D: figure 1 using ratings

ttest q1rating2 if wave==1 & inlist(q1random, 1, 2), by(q1control)
matrix mat1=(r(mu_1)-r(mu_2), r(mu_1)-r(mu_2)-1.96*r(se), r(mu_1)-r(mu_2)+1.96*r(se))
matrix rownames mat1 = a 

ttest q1rating2 if wave==2 & inlist(q1random, 1, 2), by(q1control)
matrix mat2=(r(mu_1)-r(mu_2), r(mu_1)-r(mu_2)-1.96*r(se), r(mu_1)-r(mu_2)+1.96*r(se))
matrix rownames mat2 = b 

ttest q1rating2 if inlist(q1random, 1, 2), by(q1control)
matrix mat3=(r(mu_1)-r(mu_2), r(mu_1)-r(mu_2)-1.96*r(se), r(mu_1)-r(mu_2)+1.96*r(se))
matrix rownames mat3 = c 


coefplot (matrix(mat1[.,1]), ci("mat1[.,2] mat1[.,3]") ) ///
	(matrix(mat2[.,1]), ci("mat2[.,2] mat2[.,3]") ) ///
	(matrix(mat3[.,1]), ci("mat3[.,2] mat3[.,3]") ), ///
 coeflabel(a = "Wave 1" b = "Wave 2"  c = "Wave 1 and 2") ///
mlabel format(%9.2g) mlabposition(1) xline(0, lp(dash)) graphregion(color(white)) ///
ciopts(recast(rcap)) legend(off)

**Appendix E: heterogeneous treatment effects

ttest dis if age30==0, by(q1control)
matrix mat1=(r(mu_1)-r(mu_2), r(mu_1)-r(mu_2)-1.96*r(se), r(mu_1)-r(mu_2)+1.96*r(se))*100

ttest dis if age30==1, by(q1control)
matrix mat2=(r(mu_1)-r(mu_2), r(mu_1)-r(mu_2)-1.96*r(se), r(mu_1)-r(mu_2)+1.96*r(se))*100

ttest dis if male==0, by(q1control)
matrix mat3=(r(mu_1)-r(mu_2), r(mu_1)-r(mu_2)-1.96*r(se), r(mu_1)-r(mu_2)+1.96*r(se))*100

ttest dis if male==1, by(q1control)
matrix mat4=(r(mu_1)-r(mu_2), r(mu_1)-r(mu_2)-1.96*r(se), r(mu_1)-r(mu_2)+1.96*r(se))*100

ttest dis if college==0, by(q1control)
matrix mat5=(r(mu_1)-r(mu_2), r(mu_1)-r(mu_2)-1.96*r(se), r(mu_1)-r(mu_2)+1.96*r(se))*100

ttest dis if college==1, by(q1control)
matrix mat6=(r(mu_1)-r(mu_2), r(mu_1)-r(mu_2)-1.96*r(se), r(mu_1)-r(mu_2)+1.96*r(se))*100

ttest dis if ccp==0, by(q1control)
matrix mat7=(r(mu_1)-r(mu_2), r(mu_1)-r(mu_2)-1.96*r(se), r(mu_1)-r(mu_2)+1.96*r(se))*100

ttest dis if ccp==1, by(q1control)
matrix mat8=(r(mu_1)-r(mu_2), r(mu_1)-r(mu_2)-1.96*r(se), r(mu_1)-r(mu_2)+1.96*r(se))*100

ttest dis if international==0, by(q1control)
matrix mat9=(r(mu_1)-r(mu_2), r(mu_1)-r(mu_2)-1.96*r(se), r(mu_1)-r(mu_2)+1.96*r(se))*100

ttest dis if international==1, by(q1control)
matrix mat10=(r(mu_1)-r(mu_2), r(mu_1)-r(mu_2)-1.96*r(se), r(mu_1)-r(mu_2)+1.96*r(se))*100


matrix wave12 = mat1\mat2\mat3\mat4\mat5\mat6\mat7\mat8\mat9\mat10
matrix rownames wave12 = a b c d e f g h i j

**

coefplot (matrix(wave12[.,1]), ci("wave12[.,2] wave12[.,3]") label(Wave 1 and 2) ms(Sh)),  ///
 coeflabel(a = "Younger than 30" b = "30 and Above"  c = "Female" ///
 d = "Male" e = "No College" f = "College and Above" g = "Non-CCP Member" /// 
 h = "CCP Member" i = "Little Interest"  j = "Lots of Interest") ///
mlabel format(%9.2g) mlabposition(1) xline(0, lp(dash)) graphregion(color(white)) ///
ciopts(recast(rcap))
